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ABSTRACT 

The gravitational magnification and demagnification of Type la supernovae (SNe) modify their positions on 
the Hubble diagram, shifting the distance estimates from the underlying luminosity-distance relation. This can 
introduce a systematic uncertainty in the dark energy equation of state (EOS) estimated from SNe, although 
this systematic is expected to average away for sufficiently large data sets. Using mock SN samples over the 
redshift range < z < 1.7 we quantify the lensing bias. We find that the bias on the dark energy EOS is less 
than half a percent for large datasets (> 2,000 SNe). However, if highly magnified events (SNe deviating by 
more than 2.5a) are systematically removed from the analysis, the bias increases to ~ 0.8%. Given that the 
EOS parameters measured from such a sample have a Itr uncertainty of 10%, the systematic bias related to 
lensing in SN data out to z 1 .7 can be safely ignored in future cosmological measurements. 
Subject headings: cosmology: observations — cosmology: theory — supernova — parameter estimation — 
gravitational lensing 



1. INTRODUCTION 

Since the discovery of the accelerating expansion of the 
universe dRiess et al.ll998tlPerlmutter et al.ll999HKnop et al.l 
l200llRlesret al.ll2004 . the quest to understand the physics 
responsible for this acceleration has been one of the major 
challenges of cosmology. At present the dominant expla- 
nation entails an additional energy density to the universe 
called dark energy. The physics of dark energy is gener- 
ally described in terms of its equation of state (EOS), the ra- 
tio of its pressure to density. In some models this quantity 
can vary with redshift. While there exist a variety of probes 
to explore the nature of dark energy, one of the most com- 
pelling entails the use of type la supernovae to map the Hub- 
ble diagram, and thereby directly determine the expansion 
history of the universe. With increasing sample sizes, SN dis- 
tances can potentially provide multi ple independent estimates 
of the EOS when binned in redshift (iHuterer & Cooravll2005l: 
ISuUivan. Coorav. & Holzil2007l ISulUvan et al.ll2007h. Several 
prese nt and future SN surveys, such as SNLS (lAstier et afl 
|2006|) and the Joint Dark Energy Mission (JDEM), are aimed 
at constraining the value of the dark energy EOS to better than 
10%. 

Although SNe have been shown to be good standardiz- 
able candles, the distance estimate to a gi ven SN is de- 
graded due to gr avitationa l lensing o f its flux (iFriemanll 19971 : 
IWambsgansset al. 1997; iHolz & W ald 1991). The lensing 
becomes more prominent as we observe SNe out to higher 
redshift, with the extra dispersion induced by lensing becom- 
ing comparable to the intrinsic dispe rsion (of ~ 0.1 magni- 
tudes) at z > 1 .2 jHolz & Linded2005l) . In addition to this dis- 
persion, which leads to an increase in the error associated with 
distance estimate to each individual supernova, lensing also 
correlates distance estimates of nearby SNe on the sky, since 
the lines-of-si ght pass through correlated foreground large- 
scale structure (ICooray. Huterer. & Holzl2006HHui & Greend 
l2006l) . Although this correlation error cannot be statistically 
eliminated by increasing the number of SNe in the Hubble di- 
agram, the errors can be controlled by conducting sufficiently 



wide-area (> 5 deg ) searches for SNe (in lieu of small-area 
pencil-beam surveys). 

In addition to the statistical covariance of SN distance esti- 
mates, gravitational lensing also introduces systematic uncer- 
tainties in the Hubble diagram by introducing a non-Gaussian 
dispersion in the observed luminosities of distant SNe. Since 
lensing conserves the total number of photons, this system- 
atic bias averages away if sufficiently large numbers of SNe 
per redshift bin are observed. In this case the average flux 
of the many magni fied and demagnified SNe converges on 
the unlensed value jHolz & Lindej|2005h . Nonetheless, even 
with thousands of SNe in the total sample it is possible that 
the averaging remains insufficient, given that one may need 
to bin the Hubble diagram at very small redshift intervals to 
improve sensitivity to the EOS. Furthermore, SNe at higher 
redshifts are more likely to be significantly lensed. If "ob- 
vious" outliers to the Hubble diagram are removed from the 
sample, this introduces an important bias in cosmological pa- 
rameter determination, and can lead to systematic errors in the 
determination of the dark energy EOS. 

In this paper we quantify the bias introduced in the es- 
timation of the dark energy EOS due to weak lensing 
of supernova flux. We consider the effects due to the 
non-Gau ssian nature of the lensing magnification distri- 
butions (IWang. Holz. & Munshil |2002|) . performing Monte- 
Carlo simulations by creating mock datasets for future JDEM- 
like surveys. The paper is organized as follows: In §2.1 we 
discuss our parameterization of the dark energy EOS, §2.2 dis- 
cusses gravitational lensing, and §3 is an in-depth description 
of our methodology. We present our results in §4. 

2. BACKGROUND 

2.1. Parameter Estimation Using SNe la 

For a SN with intrinsic luminosity C, the distance modulus 
is given by 

m-M = -2.51og;,) 
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where, in the framework of FRW cosmologies, the lumi- 
nosity distance di is a function of the cosmological pa- 
rameters and the redshift. We consider the two- parameter, 
time-varying dark energy EOS ( IChevalier & Pola rski 200l 
iLinder 2003) adopted by the Dark Energy Task Force (DETF) 
jAlbrec ht et al. 2006, ). We take a flat universe described by 
the cosmological parameters Q = {/z, 57,,,, wo,Wa}, with h the 
dimensionless Hubble constant, Q,„, the dimensionless matter 
density, and wq and w,, the parameters describing the dark en- 
ergy EOS w(z) = Wo + Wfl( 1 - fl). One can then express di as 
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where H(z) is the Hubble parameter at redshift z: 
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2.2. Weak Lensing of Supernova Flux 

Light from a distant SN passes through the intervening 
large-scale structure of the universe. This causes a modifi- 
cation of the observed flux due to gravitational lensing: 



irobs.lensed/^ jl\ 
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where ^{z, n) is the lensing induced magnification at redshift 
z in the direction of the SN on the sky, n, and JT"*"* '^"<= is the 
flux that would have been observed in the absence of lensing. 
The magnification ji can be either greater than (magnified) or 
less than (demagnified) one, with /i = 1 corresponding to the 
(unlensed) pure FRW scenario. 

This magnification (or demagnification) of the observed 
flux leads to an error in the distance modulus, which can be 
expressed as 



[A(m-M)],_=-2.51og 
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where we have written /i(z,n) as ^ for brevity. 

Even in the absence of lensing the measured distance mod- 
ulus to a Type la SN suffers an intrinsic error, since the super- 
novae are not perfect standard candles. This error is typically 
taken to be a Gaussian distribution in either flux or magnitude, 
with a redshift-independent standard deviation of (Tint. On 
the other hand, as lensing intrinsically depends on the optical 
depth, which increases with redshift, the scatt er (or variance) 
due t o lensing also increases with redshift dHolz & Lindej 
l2005h . The probability distribution function (PDF) for lens- 
ing magnification, P(/_t,z), of a background source at redshift 
z, depends on both the underlying cosmology and on the na- 
ture of the foreground structure responsible for lensing. We 
make use of an analytic form for P{ijl) that was cahbra ted with 
numerical simulations (IWang. Holz. & Munshill2002l) . which 
is valid for events that are not strongly lensed (/i is less than a 
few). 

The generic gravitational lensing PDF peaks 

at a demagnified val ue, with a long tail to 

high magnification (|Wa ng. Holz. & Munshil 120021: 
ISereno. Piedipalumbo. & Sazhin , ,2002) . However, since 
the total number of photons is conserved by lensing, all lens- 
ing distributions preserve the mean: (/i) = J jiPip^dji = 1. 
To ensure that this criterion is met (to an accuracy of one 
part in a million), we re-normalize the magnification PDF 
in a sligh tly different way than what was originally sug- 
gested in IWang. Holz. & Munshil ( |2002|) . P{ii) is related. 
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Redshift 

Fig. 1. — Top panel: The Hubble diagram with a lineal' redshift scale 
showing the distance modulus with redshift of a subset of a 2000 SN mock 
dataset. Bottom panel: The residuals of the data relative to the fiducial model 
of a flat universe with Q,„ = 0.3, wq = —1, w„ = 0, and /; = 0.732. 



at a given redshift, to the probability distribution for the 
reduced convergence, tj , throu gh f(/x) = Pir])/i2\Kmm\) (see 
IWang. Holz. & Munshil (l2002h : Eq. (6)), where k„„„ is the 
minimum convergence. The free parameter in normalizing 
P(i]), and hence, P(/i), is ?7,„a.t, the maximum value for the 
reduced convergence. We determine the unique value of ri,„ajt 
which yields both / P(fi)diJ. = 1 and (/i) = j iiP{^)dii = 1 
(to better than 10"^). Since the mode is not equal to the 
mean, the distribution is manifestly non-Gaussian. The 
majority of distant supernovae are slightly demagnified, and 
the inferred luminosity distances are skewed to higher values. 
The SN redshifts, on the other hand, remain unaffected, 
since gravitational lensing is achromatic. Although lensing 
magnification is insignificant at low redshifts, at redshifts 
above one both weak and strong lensing become more 
prominent. For small SN samples the overall bias is towards 
a larger cosmological acceleration (e.g., too large a value of 
r^A for ACDM models). The lensing also adds a systematic 
bias to estimates of the dark energy EOS, shifting to a more 
negative value of the EOS (e.g., less than -1.0 for a universe 
with a cosmological constant). Along with a large fraction of 
demagnified SNe, lensing produces a small number of highly 
magnified sources. If this high-magnification tail can be well 
sampled, the average of the lensing distribution is expected 
to converge on t he true mean, and the bias can be eliminated 
( iWana,200a:,Holz & Lindeill2005h . 

3. METHODOLOGY 

We generate mock SN samples, and quantify the bias in 
the estimation of the dark energy EOS due to gravitational 
lensing. We pick as our fiducial cosmology a flat universe 
with Vl,„ = 0.3, Wo = -1, Wa = 0, and h = 0.732 with a Gaus- 
sian prior of a{h) = 0.032. We consider a range of possible 
upcoming surveys, and Monte-Carlo SN samples of varying 
sizes distributed uniformly in the redshift range 0. 1 < z < 1 .7. 
If the SN intrinsic errors are Gaussian in flux, then it makes 
sense to analyze the data via flux averaging. In this case both 
the intrinsic errors and the lensing errors will average away 
for sufficiently large SN samples. However, if the intrinsic 
errors are Gaussian in magnitude, then a magnitude analysis 
may be more appropriate, although this will lead to a bias due 
to lensing (which is only bias-free for flux analysis). In what 
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Fig. 2. — Left panel: Histograms showing tlie distribution (after 10,000 realizations) of the values obtained for wq using the magnitude analysis (§ 3.1), after 
marginalizing over Wa and h, for three different sample sizes. The un-filled histogram depicts the case of the 300 SN sample, the shaded histogram shows the 
distribution for a sample size of 2000, and the hatched histogram shows the case for 10,000 SN sample. The vertical lines at -1.0095 (dashed and dotted), -1.003 
(dashed), and -1.002 (solid) show the average W(, values for the 300 SN, 2000 SN, and 10,000 SN cases, respectively. The width of the distributions is due 
primarily to the intrinsic noise of the SNe, while the shifted mode is due to finite sampling of the gravitational lensing PDF. Right panel: Histograms showing 
the distribution (after 10,000 reahzations) of the values obtained for wo using the flux-averaging technique (§ 3.2), after marginalizing over vtVi and ft, for three 
different sample sizes. The un-filled histogram depicts the case of the 300 SN sample, the shaded histogram shows the distiibution for a sample size of 2000, and 
the hatched histogram shows the case for 10,000 SN sample. The vertical lines at -1.007 (dashed and dotted), -1.0025 (dashed), and -1.001 (solid) show the vv'o 
values that were obtained after averaging over 10,000 realizations for the 300, 2000, and 10,000 sample sizes, respectively. 



follows we perform both analyses. 

3.1. Magnitude Analysis 

For each SN at a redshift z we draw a value of /i at ran- 
dom from the re-normali zed lensing PDF, P{pL,z), obtained 
from lWang. Holz. & Mun shi ( 2002), and use Eq. Q to evalu- 
ate [A(TO-M)]iensinj,. To model the intrinsic scatter we draw a 
number at random from a Gaussian distribution of zero mean 
and st andard deviation given by dint = 0.1 mag dAstier et alj 
(l2006l) find a scatter at the level of 0.13 mag). We combine 
the intrinsic and lensing noise with the "true" underlying dis- 
tance modulus to get the observed distance modulus: 

{m-Mf''\z) = {m-Mf\z)+ [Mm-M)\,,+[A{m-M)\,,,,,^ . 

(6) 

By repeating the above procedure for each of the SNe in a 
sample, we generate a mock Hubble diagram. We create 
mock surveys with different numbers of SNe: = 300, 2000, 
& 10,000; and then for each fixed number of SNe we gen- 
erate at least 10,000 independent mock samples to properly 
sample the distributions. In Figure[T]we show an example of 
one such dataset with a SNe sample size of 2000. The upper 
panel shows the Hubble diagram and the lower panel shows 
the residuals of the distance moduli relative to that of the fidu- 
cial cosmological model. 

For each mock Hubble diagram we compute the likelihood 
of a parameter set p by evaluating the x'-statistic: 



x'(P)^ 
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where <j{z\) is the error bar for the distance modulus of the ;th 
supernova, taken to be 0. 1 mag throughout. The projected bias 
in the estimation of wo can now be computed by marginalizing 
over Wa and h (keeping $7,,, fixed). 

3.2. Flux-Averaging Analysis 

Imiigl (l2000l) argues that averaging the flux, as opposed to 
the magnitudes, of observed SNe naturally removes the lens- 
ing bias, since the mean of the lensing distributions are equal 



to one in flux, but not in magnitude. We thus also analyze the 
mock data sets in flux. The flux averaging is done such that 
for each supernova we first use Eq. (|2]i to calculate the value 
of diiz), and then evaluate the fiducial flux using 



C 



4-7r[dl\z)?' 



(8) 



where we can take any a priori fixed value of C. We take 
the lensing PDF, P(fJ,), and convolve it with a Gaussian dis- 
tribution having zero mean and dispersion 0.1 (this corre- 
sponds to ~ 5% error in distance estimates). We then ran- 
domly draw from the convolved distribution, and multiply this 
number by !F^'^{z) to get the observed flux for a given super- 
nova. We repeat this process for each of the SNe to obtain 
a mock data set. We then follow the flux averaging recipe of 
^ang & Mukherjee (2004) and perform the likelihood anal- 
ysis to quantify the bias. As before we perform the above 
simulation a large (~ 10,000) number of times to get a distri- 
bution for the bias in parameter estimation. 

We further assume that our SNe samples are complete in 
the sense that there is no Malmquist bias or any bias due to 
detection effects over the redshift range considered. This is a 
reasonable assumption since we consider SNe out to a redshift 
of 1.7 and such a complete catalog is expected from the near- 
infrared imaging capabilities of the JDEM program. 

4. RESULTS AND DISCUSSION 

We first present our results for the magnitude case, as de- 
scribed in ^3.11 The left panel of Figure |2] shows histograms 
of the best-fit values of wq from the likelihood analysis, af- 
ter marginalizing over Wg and li. The empty histogram, which 
peaks at -1.009 (marked with a vertical dot-dashed line), is for 
the model with 300 SNe. The shaded histogram, represent- 
ing the 2,000 SN case, peaks at -1.003 (vertical dashed line), 
while the hatched histogram representing the 10,000 SN case 
has its peak at -1.002 (vertical solid line). These distributions 
have Icr widths of 0.016, 0.006, and 0.003, respectively. This 
scatter is primarily due to the intrinsic uncertainty associated 
with absolute calibration, and is not dominated by lensing. 
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Fig. 3. — Histograms showing tlie distribution of the values obtained for 
»'o after marginalizing over w„ and h for the 2,000 SN case (using the flux- 
averaging technique). The shaded histogram assumes that the full sample of 
2,000 SNe is used for parameter estimation. The hatched histogram shows the 
shift when outliers (SNe that are shifted above or below the Hubble diagram 
by more than 25% on either side) are removed from the sample. The bias in 
the distribution is due to the removal of highly-magnified lensing events from 
the sample. 

Without the inclusion of lensing, however, the distributions 
peak at exactly -1, and show no bias. The shifted mode gives 
us a rough idea of the bias to be expected, on average, due 
to lensing. We find that 68% of the time a random sample of 
300 SNe will have an estimated value for wq within 3% of its 
fiducial value, and this drops to 0.5% when a sample size of 
10,000 SNe is considered. 

The right panel of Figure|2]shows the same distributions as 
the left panel, but this time using the flux-averaging technique 
instead of averaging over magnitudes. The empty, shaded, 
and hatched histograms peak at -1.007, -1.003, and -1.001, 
respectively, showing the mean bias for the 300, 2,000 and 
10,000 SN cases (marked with dot-dashed, dashed, and solid 
vertical lines). With flux-averaging, we expect that 68% of 
the time a random sample of 300 SNe will yield a value of 
Wo within 2.5% of the fiducial value, and within 0.5% for a 
sample of 10,000 SNe. 

The 1(7 parameter uncertainty on wq ranges from the 20% 
level (for 300 SNe) to less than 5% (for 10,000 SNe), dwarf- 
ing the bias due to lensing. Thus, we need not be concerned 
about lensing degradation of dark energy parameter estima- 
tion for future JDEM-like surveys. We note, however, that 
our estimated bias on the EOS is larger than the lensing bia s 
of w < 0.001 quoted in Table 7 of IWood-Vasev et al. I ( 120071) . 
This is not surprising, given th eir use of the simple G aussian 
approximation to lensing from iHolz & Lindeil (|2005|) . which 
is less effective for low statistics. Nonetheless, we agree with 
their conclusion that len sing is negligible. A simil ar conclu- 
sion was also reached bv lMartel & Premadi I (|2007[) wh o used 
a compilation of 230 Tvpe la SNe (iTonrv et al. II2003I) in the 
redshift range < z < 1.8 to show that the lensing errors are 
small compared to the intrinsic SNe errors. 

We now discuss the bias which arises if anomalous SNe 
are removed from the sample. Gravitational lensing causes 



some SNe to be highly magnified, and it is conceivable that 
these "obvious" outliers are subsequently removed from the 
analysis. In this case the mean of the sample will be shifted 
away from the true underlying Hubble diagram, and a bias 
will be introduced in the best-fit parameters. To quantify this 
effect, we remove SNe which deviate from the expected mean 
luminosity-distance relation in the Hubble diagram by more 
than 25% (corresponding roughly to a 2.5ct outlier). The 
SN scatter is a result of the convolution of the intrinsic er- 
ror (Gaussian in flux of width 0.1) and the lensing PDF, and 
the outlier cutoff leads to a removal of ~ 50 SNe out of the 
2,000 SNe. These outliers are preferentially magnified, due 
to the strong lensing tail of the magnification distributions. 
The demagnification tail is cut off by the empty-beam lens- 
ing limit, and therefore isn't as prominent. The hatched his- 
togram in Figure |3] shows the distribution when events with 
convolved error greater than 2.5a are removed. The vertical 
dot-dashed line at -1.0075 shows the average value of wo ob- 
tained in this case, representing a bias in the estimate of wo 
roughly three times larger than when the full 2,000 SNe are 
analyzed (shown by shaded histogram). This bias is a result 
of cutting off the high magnification tail of the distribution, 
and thus shifting the data towards a net dimming of observed 
SNe, leading to a more negative value of wq. 

We also apply a cutoff at 3(7, in addition to the 2.5(7 dis- 
cussed above. This results in a removal of ^ 20 SNe on aver- 
age, for each 2,000 SN sample, and leads to a bias of 0.6%. 
Any arbitrary cut on the (non-Gaussian) convolved (lensing 
-I- intrinsic) sample leads to a net bias in the distance relation, 
and even for large outliers and large SN samples, this can lead 
to percent-level bias in the best-fit values for wq- 

To summarize, we have quantified the effect of weak grav- 
itational lensing on the estimation of dark energy EOS from 
type la supernova observations. With generated mock sam- 
ples of 2,000 SNe distributed uniformly in redshift up to 
z~1.7 (as expected in future surveys like JDEM), we have 
shown that the bias in parameter estimation due to lensing 
is less than 1% (which is well within the la uncertainty ex- 
pected for these missions). Analyzing the data in flux or mag- 
nitude does not alter this result. If lensed supernovae that are 
highly magnified (such that the convolved error is more than 
25% from the underlying Hubble diagram) are systematically 
removed from the sample, we find that the bias increases by a 
factor of almost three. Thus, so long as all observed SNe are 
used in the Hubble diagram, including ones that are highly 
magnified, the bias due to lensing in the estimate of the dark 
energy EOS will be significantly less than the la uncertainty. 
Even for a post-ZDZiM program with 10,000 SNe, lensing bias 
can be safely ignored. 
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